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We present an ab initio determination of tiie sliear viscosity rj of the unitary Fermi gas, based 
on finite temperature quantum Monte Carlo calculations and tlie Kubo linear-response formalism. 
We determine the temperature dependence of the shear viscosity-to-entropy density ratio rj/s. The 
minimum oi rj/s appears to be located above the critical temperature for the superfluid-to-normal 
phase transition with the most probable value being (jj/s)^;^ ~ 0.2 h/kg, which is close the Kovtun- 
Son-Starinets universal value h/{4:TTkg). 
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The unitary Fermi gas (UFG) represents a dilute but 
strongly correlated system, where the s-wave scattering 
between fermions saturates the unitarity bound for the 
cross section cr(fc) < An/k'^ [k being the relative wave 
vector of colliding particles). The system is, therefore, 
characterized by the absence of intrinsic scales, making 
it universal, i.e., independent of the details of the interac- 
tion. On the other hand, the effects of interaction have to 
be treated nonperturbatively because of the lack of any 
small parameter. The extraordinary progress in experi- 
mental methods over the last decade has brought about 
the physical realization of such a sj;stem in the form of an 
ultracold gas of fermionic atoms [1]. As a consequence, 
the UFG has provided a new paradigm for many strongly 
interacting Fermi systems, attracting attention of theo- 
retical physicists in various areas, including string theory, 
the quark-gluon plasma, neutron stars, nuclei, and to a 
certain extent high-Tc superconductivity 0]. 

Over the last few years, an impressive effort has been 
underway, both experimentally and theoretically, to es- 
tablish the physical properties of the UFG and reveal 
its strongly correlated nature. One of the most promi- 
nent manifestations of such strong correlations is the 
observation of nearly ideal hydrodynamic behavior [3- 
[^. Studies of the transport properties of these systems 
are largely inspired by a conjecture formulated by Kov- 
tun. Son, and Starinets (KSS) of the existence of a lower 
bound 77/s ^ h/{4TTkg) on the ratio of the shear viscos- 
ity rj to the entropy density s for any system [(||. As the 
bound is saturated for the case of strongly coupled Af = 4 
supersymmetric Yang-Mills theory, it is expected that 
strongly correlated quantum systems are close to this 
bound. Indeed, very different physical systems known 
to be strongly interacting appear to be very close to the 
KSS bound: i) the quark-gluon plasma created in heavy 
ion collisions at the RHIC obey rj/s ^ OAh/kg, ii) ul- 
tracold atomic gases at unitarity display rj/s ^ O.bh/kg, 
see and references therein for an extensive overview. 
It has also been predicted that low-energy electrons in 
graphene monolayers are characterized by a low value of 



77/s, of the same order as that of the quark-gluon plasma 
and ultracold atomic gases [8j]. 

In general, viscous (nonsuperfluid) hydrodynamics is 
characterized by two viscosity coefficients: the shear vis- 
cosity 77 and the bulk viscosity C. Contrary to the quark- 
gluon plasma, where the bulk viscosity is nonzero and 
can be a significant source of dissipation (especially near 
a phase transition), the bulk viscosity of the UFG van- 
ishes as a result of scale invariance [9l4ll|. The UFG is, 
therefore, an excellent candidate for a perfect fluid, de- 
fined as the one with the lowest transport coefficients rj 
and C allowed by quantum mechanics. 

A large class of theoretical methods has been used to 
determine the transport coefficients of the UFG for ho- 
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mogeneous and trapped systems 
initio calculation of the shear viscosity of the UFG is pre- 
sented within the framework of the Path Integral Monte 
Carlo (PIMC) approach [21], which has been success fully 



used to compute other properties of the UFG j22H24|. 
Contrary to the previous ab initio calculations with the 
"quenched" approximatioii, in which the fermion deter- 
minant is set to unity [2I, [Si , we compute the viscosity 
for a system with dynamical fermions. The fact that 
such a fully dynamical calculation is at all possible is not 
a priori obvious and should be regarded as one of our 
most important results. While statistical errors are ex- 
plicitly under control, we provide only a limited assess- 
ment of systematic effects (finite density and volume). 
From our results it is clear that those effects can be con- 
trolled. While we focus our study on the shear viscosity, 
we have preliminary indications that the bulk viscosity 
vanishes at all temperatures, in agreement with the scale 
invariance arguments mentioned above. However, we de- 
fer more careful determinations of both viscosities as well 
as better control of systematic errors to future work. 

Transport coefficients can be theoretically determined 
using hnear response theory via the Kubo relations 
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In order to apply such relations within the 
framework of PIMC calculations, we followed the method 
based on the stress-tensor correlators [25, M,^- Within 
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this approach, the frequency-dependent shear viscosity is 
given by (in units such that h = kg = m = 1) 
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(1) 



while the static viscosity is defined in the limit of zero 
frequency: r/ = lim;^_^o+ vi^)- The spectral density 
Pijkiil^^) related to the imaginary-time (Euclidean) 
stress-tensor correlator G^j {q, r) by inversion of the re- 
lation 

r f , r ( ^ coshKr-/3/2)] ^ 
G.,,uM,r)^j^ p,,,,,{q,u) du;, (2) 

where (3 = 1/T is the inverse temperature. In turn, the 
stress-tensor correlator has the form 

G„- fe, (q, t)= j d^re-^"-^ {% (r , r)II,, (0, 0)) , (3) 

where the average is performed over the grand canon- 
ical ensemble, 0(r) = ^r{H~f,N)Q^-r{H-,,N) ^ fj jg ^-^^ 

Hamiltonian of the system, fi is the chemical potential, 
and N is the particle number operator. The stress-tensor 
operator 11^^ (r) is defined via the operator version of 
the Euler equation (summation over doubled index is as- 
sumed) : 



iMr),H]^din^i{r), 



(4) 



where jj^ is the current operator. Since the current 
operator commutes neither with the kinetic-energy nor 
with the potential-energy parts of the Hamiltonian, it 
is convenient to split the stress tensor into two parts: 



n 



kl 



n 



(T) 



kl + ^^ki^ ■ The kinetic-energy part IIj,; 
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is well 



established and is the only contribution to the shear vis- 
cosity for a zero-range potential (see for example [lit). 
The potential-energy part H^]'^ is more complicated, as 
defining the diagonal of the stress tensor is not trivial 
due to scale invariance, which is violated in our lattice 
calculations. Nevertheless, if we proceed with the stress 
tensor which on the lattice does not respect the sum rule 
J (Prll^^(r) = 2H imposed by the scale invariance [lo| . 
we obtain results consistent with C = 0- This matter is 
under further investigation. 

Using the PIMC method, the stress-tensor correla- 
tor ([3]) was evaluated at g = for 51 points in imaginary 
time T, uniformly distributed in the interval [0,/3] on a 
spatial lattice of 8'^ points. Increasing the number of r 
points did not affect the final results. A statistical ensem- 
ble of 5000 uncorrelated samples was generated at each 
temperature, thus reducing the statistical errors to a few 
percent (depending on the temperature and value of r). 
To estimate the size of discretization errors, exploratory 
calculations on a 10'^ lattice were performed. All the cal- 
culations presented here were performed with an average 
particle number density n = N/V w 0.09. The system- 
atic errors associated with the stress-tensor correlator. 



related to finite volume effects as well as effective-range 
corrections, are likely ~ 10 — 15% 21, 2^. For a more 
detailed discussion, see Ref. [30| . 

To determine rj, one has to solve Eq. ([1} numerically, 
which is an ill-posed inversion problem, as there exist 
an infinite number of solutions that reproduce the cor- 
relator within its error bars. Therefore, estimating the 
shear viscosity requires additional information. Besides 
the non-negativity of the viscosity ri(uj) ^ 0, the sum 
rule and the asymptotic tail behavior (see |ll| with sub- 



sequent corrections [15|, |3l|) have been used as a priori 



information. In the unitary limit these conditions read 
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where C is Tan contact density [32] and e is the energy 
density. The energy density is obtained directly from 
PIMC calculations, while the contact density is taken 
from Ref. [33]. Based on the results for the nonintcr- 
acting Fermi gas, where 77pQ(aj) oc S (uj), and those ob- 
tained within the T-matrix approach |15| or kinetic the- 
ory [l^, the shear viscosity ri{uj) is expected to be a 
continuous function with Gaussian-like structure at low 
frequencies, smoothly evolving into the asymptotic tail 

C 

behavior r]{Lu co) ±2 — ; Moreover, we assume 



that there is no sharp structure in the spectral density 
in low frequency limit (associated, for example, with well 
defined quasiparticles) , which could be overlooked dur- 
ing the inversion process. We used these assumptions to 
construct the model used in the inversion procedure. 

To perform the inversion we applied a methodology 
based on two complementary methods: singular value 
decomposition (SVD) and maximum entropy method 
(MEM), both described in Ref. (SJ]. Since these methods 
are based on completely different approaches, a solution 
that is in agreement simultaneously with both of them is 
regarded as the most favorable scenario. In order to esti- 
mate the stability of the combined methods with respect 
to the algorithm parameters, the "bootstrap" strategy 
was applied. Namely, about 200 reconstructions were 
performed, with randomly generated initial parameters 
(within some reasonably chosen interval). The collected 
set of samples was subsequently used to evaluate the av- 
erage value of the shear viscosity and the standard devi- 
ation (see [sOl for details). 

In Fig.[l] the dimensionless static shear viscosity rj/n is 
shown as a function of T/sp, where Sp — (37r^n)^/'^/2m 
is the Fermi energy of the noninteracting gas. The shear 
viscosity monotonically decreases with decreasing tem- 
perature. No drastic suppression of the viscosity below 
the critical temperature of the superfluid-normal phase 



transition ~ 0.15£^ is observed. However, note that 
below Tj, the coefficient rj describes the viscosity of the 
normal fluid component only. The results on 8^ and 10'^ 
lattices exhibit satisfactory agreement. Surprisingly, our 
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FIG. 1: (Color online) The dimensionless static shear viscos- 
ity ri/n as a function of T /ep for an 8'' lattice (red) squares 
and 10"^ lattice solid (blue) circles. The error bars only rep- 
resent the stability of the combined (SVD and MEM) inver- 
sion procedure with respect to changes in the algorithm pa- 
rameters. The (green) line depicts the prediction of kinetic 
theory [13]. For comparison, recent results of the T- matrix 
theory produced by Enss et al., are plotted as open (purple) 
circles llal. 



results approach the predictions of kinetic theory already 
at r > O.Se^ ji2l- Note that the PIMC resuhs are sig- 
nificantly below all known results in the vicinity of T^. 

In Fig.[2l the value of the entropy obtained from PIMC 
calculations is shown (extracted as in Ref . [2l| ) , together 
with the results extracted from the recent high-precision 
MIT measurement [s^. For temperatures T > 0.25£^, 
both lattices reproduce experimental data reasonably 
well. At low temperatures T < 0.25e^ the 8'^-lattice 
results deviate from the measurements, producing sys- 
tematically lower values. On the other hand, the 10'^- 
lattice results reproduce correctly the temperature de- 
pendence of the entropy, yet slightly overestimating the 
experimental values. These discrepancies are attributed 
to systematic errors that are known to be present at low 
temperatures even for larger lattices [IS]- Consequently, 
we expect the ratio rj/s to be significantly affected by 
uncertainties related to the entropy at low temperatures. 

In Fig.[3]the ratio rj/s is presented as a function of tem- 
perature. The PIMC calculations reveal the existence of 
a deep and rather narrow minimum in rj/s at tempera- 
tures around 0.20 — 0.25ej^, which is above T^. Again, the 
ratio r]/s is located around the kinetic theory predictions 
already at T > O.Se^ fi^. The estimation of the 77/s- 
ratio reveals (?7/s)j„i„ ~ 0.2 as the most probable value 
for the minimum. This result is about 2.5 times higher 
than the KSS bound 77/s > l/47r « 0.08. Such a low 
value has been reported o nly for pure gluons as a result 
of lattice calculations [2^, [2g| . 

The minimum value for the ratio ('7/s)min ~ 0.2, is 



FIG. 2: (Color online) Entropy per particle as a function of 
T/sp for the 8'^ lattice in (red) squares and 10'^ lattice in 
(blue) circles. The entro py p er particle extracted from the 
recent MIT measurement [33 is plotted with (black) crosses. 



significantly lower than predictions of all current cal- 
culations, which yield a minimum ~ 0.5. However, 
these methods are in principle unreliable when applied 
to the UFG at T ~ T^, where the minimum appears. 
Moreover, the rj/s ratio calculated from PIMC simula- 
tions is also significantly lower than the experimental 
measurements [SnSj, which also give the value ~ 0.5. 
Note, however that these measurements are performed 
in trapped systems. The trap-averaged viscosity {rj/n) — 
Wh I ^(''") affect the determination of the min- 

imum value. To solve this puzzle, one should apply an 
averaging procedure to the uniform case results, using, 
e.g., local density approximation. It is well known that 
this procedure leads to a divergence due to the viola- 
tion of the hydrodynamic description at the edges of the 
cloud [sgI . To perform a reliable averaging procedure the 
coUisionless edges should be treated using kinetic theory. 
This, however, is a hard task that requires the knowledge 
of second-order transport coefficients like the relaxation 
time, which are currently poorly known. 

Since our main result for the minimal value of 77/s is 
significantly lower than other predictions as well as ex- 
perimental results, we have performed exploratory cal- 
culations to estimate the size of systematic effects. We 
have checked the stability of the inversion procedure with 
respect to the default model as well the impact of the 
nonzero value of the effective range, see [s^ for details. 
Our conservative estimation indicates that the minimal 
value of the ry/s-ratio is lower than 0.45. 

In summary, we have presented an attempt to deter- 
mine the shear viscosity of the UFG through an ab ini- 
tio PIMC approach. The minimum value of the 77/s ra- 
tio was estimated to be lower than 0.45 with the most 
probable value being {11/3)^^^ w 0.2, located around 
T « 0.20 - 0.25e^. This value is close to the KSS bound 
and suggests that the unitary Fermi gas is the best can- 
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FIG. 3: (Color online) Ratio of the shear viscosity to en- 
tropy density 77/s as a function of T/sp for an 8^ lattice (red) 
squares and 10'' lattice (blue) circles. The error bars only 
represent the stability of the combined (SVD and MEM) in- 
version procedure with respect to the change of algorithm pa- 
rameters and do not include systematic errors of the entropy 
determination. Results of the T-matrix theory are plotted by 
open (purple) circles [31 . In the high- and low-temperature 
regimes, known asymptotics are depicted: for T > 0.3ep 
the prediction of kinetic theory 112] as a green line, and for 
T < 0.2ep the contribution from phonon excitations [l^l as a 
brown line. The KSS bound appears as a dashed black line. 



didate for the perfect fluid. As our results can be signif- 
icantly affected by systematic errors, further and more 
precise investigations are called for. 
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ANALYTIC CONTINUATION & INVERSION 

The calculation of dynamic response functions, such 
as susceptibilities, viscosities and conductivities, entails 
computing real-time correlation functions. Since PIMC 
calculations are performed in imaginary time, one faces 
the problem of analytic continuation. As real-time and 
imaginary-time correlations share the same spectral den- 
sity p, this problem can be recast as an inversion problem 
in which one attempts to find p with imaginary-time cor- 
relation functions as the starting point. In the case of 
the shear viscosity this inversion problem is given by the 
equation (for zero momentum q = 0): 



G^y,^y{Ti)^ P^y^^y{uj)K{uj,T^)duJ, (1) 

Jo 



with the kernel 



KiLU,T) 



cosh (w(t- /3/2)) 
sinh M/2) 



(2) 



The correlator G^y^y{T^) is determined within PIMC 
with certain accuracy /S.G^y^y{T^) for a finite set of 
points T]^, . . . , uniformly distributed within interval 
[0, 13] (we used iV = 51 for 8^ and iV = 41 for 10^ lattice). 
For brevity we skip subscripts (xy, xy) hereafter. By def- 
inition, the spectral density vanishes at zero frequency, 
while the kernel has a pole K{{), r) — )■ oo. Consequently, 
the expression p{Q)K{Q,t) is not a well defined quantity 
for the lower limit of the integral ([T]). However, since we 
are interested in the static shear viscosity 



lim ii{uj) = lim 



■Kp{uj) 



(3) 



we can formulate the problem ([T]) as 
1 



G(t,) = - / ri{uj)K{Lo,T,)du, (4) 



with a new kernel K{uj,t) — ujK{uj,t). Note that the 
new kernel is well defined at zero frequency K{0,t) = 
2/13. Moreover, the frequency dependent shear viscos- 
ity can be now directly determined, without taking the 
limit lim(^_j.Q+ p{lj)/uj which could be difficult to realize 
numerically. 

From the symmetry of the kernel we infer a simple 
relation for the correlator G{t) — G{/3 — r), which is 
simply a result of the bosonic character of the stress 



tensor, and is fulfilled in our calculations within er- 
ror bars. Consequently, we may restrict the inversion 
to the interval t € [0,/3/2]. Moreover, we have re- 
placed G(r) 4- [G{t) + G{l3 - t)]/2 in order to obtain 
"smoother" , symmetric correlator. The problem ([1} is 
supplemented with external constraints: non-negativity 
of the shear viscosity, sum rule and asymptotic tail be- 
havior. From the known tail behavior we obtain that the 
correlator G(r) should have poles for t e {0,(3}. Indeed, 
the contribution from the tail for r = reads (in units 



such that h ^ kg ~ m = 1) 



+ 00 



G 



large 



' cLu = 



(5) 



where uj^^^^^ 1 is assumed to be large enough to ap- 
proximate the viscosity by the tail behavior and K(uj, 0) 
by unity. The above result indicates that the the correla- 
tor at small r does not carry much information about the 
low-frequency part of the shear viscosity. The most im- 
portant information about the shear viscosity at low fre- 
quencies is encoded in the region located around t = (3/2. 
Moreover, the PIMC approach does not provide us with 
a correlator G(r) that acquires extremely high values at 
the edges of domain. We attribute this anomaly to the 
systematic error related to the existence of the cut-off in 
the momentum space, which implies a cut-off in the fre- 



quencies at u!„ 



Pmax/2- Since for small r, G(r) is 



strongly affected by systematic errors and does not pro- 
vide significant information for estimation of the static 
shear viscosity, we remove from the analysis few first 
points of G(rj) (typically about 3). 

The inversion was performed using a combination of 
two complementary methods: Singular Value Decompo- 
sition (SVD) and Maximum Entropy Method (MEM). 
We have applied a variant of the MEM referred to as self- 
consistent MEM in the paper [l[. Within this method, 
the a priori model solution is not fixed (as it is in stan- 
dard MEM) but evolves simultaneously with the solution 
within the space of functions spanned by admissible mod- 
els. 

The class of models for the self-consistent MEM is de- 
fined as 



M(a;, {m, a, c, «!, aa}) = /(w, {a^, q;2}): 



G 



+ [l-/(w, {a„a^})]N{u, {m, a, c}), (6) 



— where 



N{u;, {m, a, c}) = ^=L= exp f - ^^^^^-J^ ) , (7) 



and 



fiuj,{a^,a^}) = e 



e«i" _ 1 

1 gQiCw-Qj) ' 



(8) 
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FIG. 1: (Color online) Fit of the proposed model (blue line) 
to the results produced by T-matrix theory (red and purple 
symbols) ^ for two selected temperatures. 



such that f{uj ^ 0) ^ and f{uj — > oo) — > 1. The func- 
tion /(a;, {a^, Q!2}) guarantees a smooth change of be- 
havior between Gaussian dependence and the known tail 
dependence. The set of five parameters {m, a, c, ai, 
describes the available degrees of freedom of the model 
for the self-consistent MEM. To initialize the model for 
the first iteration we fit it to the SVD solution. We 
have found that the proposed model can reasonably de- 
scribe the SVD solution for all considered cases. We have 
also checked that the solutions provided by the T-matrix 
theory [2] can be well reproduced by the used model, 
see Fig. [TJ We also applied the T-matrix results, to esti- 
mate the range of frequencies uj^^^^ above which the solu- 
tion is well reproduced by the tail behavior. 

The uncertainty related to the inversion procedure was 
estimated using the "bootstrap" strategy. The bootstrap 
sample consists of about 200 launches of the algorithm 
with randomly generated (from some reasonable inter- 
vals) parameters: i) ui^ax ~ upper limit of integra- 
tion, ii) tj-'taii " point where the universal tail behavior 
starts, iii) a - the parameter for the MEM algorithm. To 
represent the MEM solution r]{uj), we used mesh of 100 
points uniformly distributed within range ui S [0,i^„iax]- 
The static shear viscosity 77(0) is computed as the aver- 
age over the bootstrap collection, while the uncertainty 
is determined by the standard deviation. 

In Fig. [2] we present a sample of results for three tem- 
peratures. For temperatures above 0.25e^ we observe 
that 77 (cj) possesses a Gaussian-like structure at low fre- 
quencies. This structure is smeared out at low tempera- 
tures. Note also that for low temperatures both the SVD 
and MEM solutions are similar, while at larger tempera- 
tures the SVD method produces solutions of lower static 
shear viscosity 77(0) than the MEM approach. Such dis- 
crepancy is permissible, since the SVD method gives only 
the projection of the solution onto the relatively small 
subspace where the inverse problem is well defined [T] . 



0.6 
0.5 
0.4 

1; 0.3 

0.2 
0.1 


0.7 





SVD 




MEM 




initial modei 




finai modei 




taii asymptotics 









0.3 
0.2 
0.1 





1 2 3 4 5 6 7 





SVD 




MEM 




initial model 




final model 




tail asymptotics 









1 2 3 4 5 6 7 

(O/Ep 



FIG. 2: (Color online) The frequency dependent shear vis- 
cosity in units of density for three different temperatures ob- 
tained from the Quantum Monte Carlo calculations with 8^- 
lattice. The SVD solution is used to create the initial model 
for the self-consistent MEM. 



DISCUSSION OF SYSTEMATIC ERRORS 

In this section we present exploratory estimations of 
corrections arising from systematic errors. Since the re- 
sults for two different lattice sizes = 8 and 10 agree 
sufficiently well, we shall focus on the corrections at- 
tributed to the inversion procedure and effective range. 

According to our experience the methodology based 
on the combination of the SVD method and the self- 
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FIG. 3: (Color online) The frequency dependent shear vis- 
cosity in units of density obtained from the self-consistent 
MEM when initialized with exaggerated model. 



consistent MEM provides the most reliable results. How- 
ever one may ask the following question: What is the 
maximal value of the shear viscosity which is simultane- 
ously consistent with QMC data and with the external 
constraints? In order to answer this question we have 
applied the self-consistent MEM with the initial model 
predicting a significantly higher value of the viscosity 
than all known predictions for the studied temperatures, 
see Fig. [3l By construction, the initial model does not 
satisfy the sum rule. Since we start with this exaggerated 
model in each subsequent iteration the self-consistent 
MEM systematically reduces r]{Lo) and thus provides us 
with an upper bound for the viscosity. In general, we ob- 
serve that the produced solutions are in clear disagree- 
ment with the solutions provided by the SVD method. 
The estimated upper bound for rj/s reveals the value 
('7/^)min ^ 0-45 (see Fig. H]), where we used the smallest 
obtained value for the entropy density. Namely, we used 
8'^-lattice results, which clearly are affected by systematic 
errors, especially in the regime where the minimum is lo- 
cated (see Fig. 2 in the main paper), which artificially 
enhances rj/s. It should also be noted that the generated 
values for the static viscosity 77 significantly overestimate 
all known results for uniform systems for temperatures 
T > OAsp. As a consequence, it is difficult to obtain 
a smooth connection of the upper limit predictions with 
the known behavior of the viscosity at high temperatures 
ri/n w 2.n{T/epf/'^ '2]. This problem is absent for the 
most probable solutions, where the static shear viscosity 
approaches smoothly results of kinetic theory. 

To obtain stable results for the inversion procedure 
with respect to change of the algorithm parameters it 
is necessary to include external constraints, especially 
the sum rule together with the asymptotic tail behav- 
ior. Both of these constraints as ingredient include Tan's 
contact density. Since we were not able to determine the 
value of the contact within present work from the asymp- 
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FIG. 4: (Color online) The dimensionless static shear viscos- 
ity r]/n (upper panel) and the ratio of the shear viscosity to 
the entropy density r;/s (lower panel) as a function of T/sp. 
Filled (red) squares and filled (blue) circles show the results 
for the 8^ and 10^ lattices, respectively, with n ~ 0.09 as pre- 
sented in the main paper. Results for reduced densities are 
shown as open (green) triangles and open (purple) diamonds. 
The (red) dotted line shows a conservative estimate for the 
upper bound. For the 77/s ratio, the upper bound includes 
also possible corrections arising from systematic errors in the 
estimation of the entropy. The dashed (black) line shows the 
KSS bound. 



totics of the momentum distribution (which requires very 
low densities), we used values obtained in earlier work [3] . 
This may introduce additional systematic error into our 
calculations, but we have checked that variation of the 
contact by ±10% generates change of the shear viscosity 
within plotted error bars. 

Recent results 043 suggest that we should expect sig- 
nificant modifications arising from effective range correc- 
tions as the product r^gkp ~ 0.6 remains non-negligible. 
To check the impact of these corrections we performed 
exploratory calculations with reduced density. For the 8"^ 
lattice we reduced the density to rt ~ 0.07 (we decided 
not to reduce it further to avoid shell effects) and for the 



10 lattice to n ~ 0.04 which corresponds to r^^kp « 0.5 
and 0.4, respectively. Performing calculations at signif- 
icantly lower densities, while keeping T/ep sufhciently 
low, is an extremely time-consuming task. We therefore 
decided to perform calculations only for a few selected 
temperatures. Our limited number of tests indicate that 
decreasing finite-range effects yields slightly lower vis- 
cosities for temperatures T/sp > 0.3, see Fig. |31 In the 
regime of temperatures where the minimum of 77/s-ratio 
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is located, all the results agree within the inversion er- 
ror bars. Hence, we conclude that the estimated upper 
bound for the shear viscosity is rather conservative. All 
the theories that produce viscosities above our estima- 
tions of the upper bound are in clear disagreement with 
our results. 
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